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ABSTRACT 

We examine some implications of inertial range and dissipation range 
correlation and spectral analyses extracted from 33 intervals of Wind magnetic 
field data. When field polarity and signatures of cross helicity and magnetic 
helicity are examined, most of the datasets suggest some role of cyclotron 
resonant dissipative processes involving thermal protons. We postulate that an 
active spectral cascade into the dissipation range is balanced by a combination of 
cyclotron-resonant and non-cyclotron-resonant kinetic dissipation mechanisms, 
of which only the former induces a magnetic helicity signature. A rate balance 
theory, constrained by the data, suggests that the ratio of the two mechanisms 
is of order unity. While highly simplified, this approach appears to account for 
several observed features, and explains why complete cyclotron absorption, and 
the corresponding pure magnetic helicity signature, is usually not observed. 

Subject headings: MHD — turbulence 
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1. Introduction 

The solar wind plasma displays many characteristics that can be reasonably well 
described by a magnetohydrodynamic (MHD) fluid model (Tu & Marsch 1995; Burlaga 
1995), including features that appear to be related to fluid turbulence (Coleman 1968) and 
MHD wave activity (Belcher & Davis 1971). Within the context of a simple nonlinear MHD 
theory, one expects that a key feature is spectral cascade of energy from larger, energy- 
containing scales through an inertial range and ultimately into a dissipative range (von 
Karman and Howarth 1938; Batchelor 1970; Martinez et al. 1997). An MHD description of 
the solar wind or other collisionless plasmas in astrophysics is a drastic oversimplification, 
and it is therefore significant that solar wind observations support the general picture of a 
turbulent MHD cascade from large to small scales (Jokipii 1973; Matthaeus &; Goldstein 
1982; Goldstein, Roberts & Matthaeus 1995). Indeed, unless turbulent transfer and decay 
are invoked, it is difficult to explain the proton temperatures at 1 AU (Coleman 1968), the 
solar wind's general non-adiabatic temperature profile (Freeman 1988; Richardson et al. 
1995) and the radial variation of the fluctuation levels (Zank, Matthaeus & Smith 1996). 

There is an appealing simplicity in an explanation of these features in terms of a 
cascade controlled by a larger scale energy-containing fluctuations (Matthaeus et al. 1994; 
Zank et al. 1996) and mediated by a self-similar inertial range (Tu, Pu & Wei 1984; Verma, 
Roberts & Goldstein 1995). This model is in direct analogy with hydrodynamic turbulence 
(Batchelor 1970). Nevertheless, there is an essential and theoretically challenging piece 
that is missing in this picture. The cascade must terminate at small scales, possibly in a 
dissipation range in which processes occur that convert MHD fluctuation energy into plasma 
thermal energy. The possible involvement of ion cyclotron activity in the observed onset 
of steepening of solar wind magnetic field spectra (at pa 1 Hz) (Behannon 1976; Denskat, 
Beinroth & Neubauer 1983) has been discussed for some time. Generally, discussion of the 
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collisionless damping of interplanetary fluctuations has concentrated on Landau damping 
(Barnes 1966, 1979). It is only recently (Goldstein, Roberts & Fitch 1994; Leamon et al. 
1998) that attention has begun to focus on a broader framework for explaining dissipation 
processes. A theoretical perspective that invokes kinetic theory to convert fluid scale energy 
to heat is needed, taking into account spectral transfer that continually resupplies the 
dissipation range through broad-band nonlinear couplings. This paper provides a simple 
description of this process based upon the assumption of a steady cascade, with the goal of 
explaining recently described features of the dissipation. 



It is useful to adopt a leading order description based upon incompressible turbulence, 
in view of the low level of interplanetary density fluctuations (Roberts et al. 1987), the 
observed density spectrum (Montgomery, Brown & Matthaeus 1987), and the low average 
turbulent Mach number (Matthaeus, Goldstein & Roberts 1990). This perspective is 



spectrum. Neglecting small internal energy fluctuations, the turbulent energy per unit 
mass, E, consists of contributions from the turbulent (ion) velocity v and the fluctuating 
component of the magnetic field b, scaled to Alfven units. For an appropriately defined 
ensemble average (•••), the contribution to the energy from velocity fluctuations E v and 
from magnetic fluctuations is 



2. MHD turbulence parameters 



also consistent with the persistence of the k 5//3 signature of the Kolmogoroff cascade 




(1) 



In its idealized definition, the turbulent energy includes contributions from all 
wavenumbers and frequencies. However, in some circumstances one might consider only 
contributions from certain scales, so that, for example, the spectral decomposition of 
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magnetic energy, E b = J d 3 kE b (k) might include only a certain range of wavenumbers. One 
might choose to look at the energy in a finite band of wavenumbers or frequencies, for 
example, when the physics of the inertial- or dissipation range is discussed. 

Apart from energy, other quantities of importance for MHD turbulence are the 
magnetic helicity H m = (b • a), where b = V x a, the cross helicity H c — (v • b), and 
the respective spectral decompositions (Matthaeus & Goldstein 1982). The amounts of 
cross helicity and magnetic helicity relative to the energy are conveniently measured by the 
following dimensionless parameters. The normalized cross helicity 



E + -E„ 



(2) 



c E + + E_ ' 

is defined in terms of the Elsasser energies E± = (|v ± b| 2 ) (Marsch & Mangeney 1987), 
and lies between —1 and +1. Normalized magnetic helicity 

° m = F I F ' ( 3 ) 

is written here in terms of El, the magnetic energy in left-handed (positive helicity) 
spatial structures, and Er, the magnetic energy in right-handed (negative helicity) spatial 
structures. Note that E b = E L + E R . We use the following sense of circular polarization: 
right-handed means a sense of rotation from the x direction towards the y direction as 
one samples in the positive z direction for a right-handed (x, y, z) coordinate system. In 
terms of the integrated magnetic helicity spectrum, E b L = \ (E b + / d 3 k\k\H m (k)) and 
E b R = \(E b -jd 3 k\k\H m (K)). 

The magnetic helicity is important in the present context because spatial handedness is 
related to resonance conditions with charged particles. Cross helicity relates to the direction 
of propagation of large amplitude Alfven waves with respect to a uniform or slowly varying 
background magnetic field B (Belcher & Davis 1971; Matthaeus & Goldstein 1982). Both 
together determine the polarization of the waves in the plasma frame (Smith et al. 1984). 
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3. Observations 

In a recent study, Leamon et al. (1998) described properties of the interplanetary- 
dissipation range at 1 AU. Their analysis included spectra and other parameters computed 
for 33 intervals of high time resolution (up to 22 vectors/s) Wind magnetic field data, along 
with plasma data at a much lower sampling rate (either 46 or 92 seconds per measurement). 
In this analysis, the magnetic field data provides information about Eb, El and Er in 
the dissipation range and in the inertial range. For the samples in the Leamon et al. 
study, the inertial and dissipation ranges were distinguished according to spectral slope. 
The average inertial range spectral index corresponded to a one-dimensional spectral law 
in good agreement with the Kolmogoroff value, Eb(k r ) ~ fc~ L67 , for radial wavenumber 
k r = 2nf /Vswi with solar wind speed V$w an d / the spacecraft-frame frequency. The 
dissipation range spectra were steeper, averaging Eb{k T ) ~ k~ 3m , with a breakpoint between 
the two ranges at an average frequency of about 0.5 Hz. 

Leamon et al. noted that most of the intervals they examined showed a signature in 
the magnetic helicity at dissipation range frequencies, as had been reported previously 
by Goldstein et al. (1994). In contrast, typical inertial range magnetic helicity spectra 
oscillate randomly as a function of frequency (Matthaeus & Goldstein 1982). Leamon et 
al. found that as much as 90% of the energy to be carried by waves propagating at highly 
oblique angles or quasi-two-dimensional turbulence rather than parallel-propagating Alfven 
waves. Nevertheless, in almost all of the intervals examined, the dissipation range H m were 
consistent with absorption of outward-propagating Alfven waves by resonant coupling to 
thermal protons. 

EDITOR: PLACE FIGURE HERE. 
Here we examine in greater detail the data underlying the latter conclusion. In Fig. 1 we 
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show the normalized cross helicity a c computed from inertial range data, plotted versus 
the normalized magnetic helicity a m in the dissipation range, for the 33 data intervals 
previously analyzed. H c can be computed only in the inertial range due to limited sampling 
rates for plasma data; we use the inertial range H c as a proxy for the same quantity that 
is unmeasurable in the dissipation range. In effect, we are assuming that the direction 
of propagation of fluctuations in the dissipation range is the same as the direction of 
propagation of fluctuations in the inertial range. 

It is apparent from the data in Fig. 1 that most intervals for which the mean magnetic 
field is outwardly directed have a m > and u c < 0. On the other hand, inwards directed 
B is associated with a m < and a c > 0. This implies a predominance of outward 
propagating waves. One can readily see that this is consistent with cyclotron-resonant 
absorption of outward-propagating fluctuations by thermal protons, as follows. A proton 
moving outward along the magnetic field executes a left-handed helical trajectory. Waves 
propagating outward at the Alfven speed will overtake most thermal particles (at (3 1) 
and therefore, on average, the thermal protons will be in resonance with such waves that 
have a right-handed spatial handedness (negative H m ). If the energy of these waves is 
assumed to be damped by the resonant protons, the energy that remains will preferentially 
reside in the undamped fluctuations, which have a left-handed structure and positive H m 
(see, for example, Moffatt (1978)). Consequently, outward B should be associated with 
a c < (outward propagation) and a m > 0. Reversing the direction of B but maintaining 
the assumption of outward propagating waves (now a c > 0) produces the conclusion that 
a m < in the dissipation range by the same argument. 
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4. Cascade and dissipation 

The above argument explains the clustering of the observational points in the upper 
left and lower right quadrants. However, there are questions that arise. First, if kinetic 
processes are assumed to be very rapid, why is the signature in the magnetic helicity not 
pure (±1) as one would expect for complete cyclotron absorption? Second, how is the 
above argument modified if instead of pure cyclotron-resonant absorption processes, there 
is also a contribution due to Landau resonance or nonresonant absorption? Finally, since 
the observed cross helicities are not "pure," what is the effect of relaxing the assumption of 
purely outward traveling Alfven waves? 

It turns out that these questions can be addressed, in at least a preliminary fashion, 
by postulating a cascade and associated dissipation processes that are described by a set of 
energy balance equations, as follows: 

dE L S 



dt 2 
dE R S 



70-El - l r E L 



The energies in left- and right-handed spatial structures are respectively designated as 
El and Er following our earlier discussion (in this case the integration over the spectrum 
now includes, by assumption, only the dissipation range). The rate of supply of energy 
(per unit mass) transferred into the dissipation range from the inertial range is designated 
by S. This supply rate is equally apportioned to L and R fluctuations since inertial range 
H m is random. We assume that the only external contribution to dE L> R/dt is due to the 
cascade term S, and that in the dissipation range there is no exchange between El and Er, 
or exchange between kinetic and magnetic energies. The quantity 7 r represents a decay 
rate due to cyclotron-resonant absorption by thermal protons, and it appears only in the 
L equation under the assumption that fluctuations are outward propagating and B is 
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inward. (This would also occur for inward propagation and outward B .) The remaining 
damping term, 70, appears in both L and R equations and represents decay processes 
that produce no signature in the magnetic helicity. Included in 70 are contributions from 
Landau damping and other mechanisms that do not involve cyclotron resonance, as well as 
mechanisms that are fully nonresonant. 



5. Cyclotron-resonant and other forms of dissipation 

We can now proceed to estimate a typical relative strength of cyclotron-resonant and 
non-cyclotron resonant processes. Supposing the cascade is steady, so dE^/dt = 0, and 
we may equate the right hand sides of Eqs. @. From the data, we take a typical value 
of magnetic helicity to be a m ~ —1/3. This corresponds in Eq. (|3|) to E R = 2E L in the 
dissipation range. Then for consistency with Eqs. (f|) we must have 70 ~ 7 r , indicating that 
cyclotron and non-cyclotron absorption mechanisms are approximately of equal strength. 

Since observed values of H c are not pure, the above argument should be refined to 
account for a distribution of propagation directions relative to the slower thermal protons. 
Assume, then, that there is a probability P(L) that fluctuations are propagating outward, 
which produces a resonance between left-handed structures and thermal protons, and 
implies the appearance of 7 r in the El equation. Assigning the probability of inward 
propagation to be P{R) = 1 — P{L) implies that resonance between right-handed structures 
and thermal protons is weighted accordingly. Therefore, the cascade balance equations 
become: 

dE L S 



dt ~ 2 



lo E L - P[L)lrE L 



^ = \-l*E R -P{R) lr E R (5) 
According to the Elsasser representation, fluctuations with energy E- tend to propagate 
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along the mean field B while fluctuations having energy E + tend to propagate antiparallel 
to B . We assume for simplicity that the probability that, at any location in the plasma, 
a typical thermal proton will "see" outward propagation is proportional to the average 
outward propagating energy. Thus, 

and therefore P(R) = (1 - a c )/2. 

With this interpretation, we can make use of the data in Fig. 1 to constrain our model 
and arrive at further insights about the dissipation processes. We invoke the steady form 
of Eqs. (H) along with the definitions Eqs. @, @ and @, and assume that 70 and 7 r are 
independent of a c , a m and other plasma turbulence parameters. Eliminating El and Er, 
we conclude that 

ac = - U + <r m . (7) 

The best-fit line forced through the origin is a c = — 1.90a m , while the best-fit straight 
line through the data is o c = — 1.80er m + 0.10. Considering either 32 or 31 degrees of freedom 
accordingly, the reduced chi-squared values of the two fits are Xr = 1-78 and 1.55. Putting 
a c = — 1.90er m in Eq. (|7|) implies that 7 r = 2.227 . The other important consequence of 
Eq. (0) is that only when 70 = do pure Alfven waves lead to purely helical states. 



6. Interpretation & Discussion 

This preliminary attempt to understand the observed interplanetary dissipation 
range spectra, while clearly oversimplified, appears to contain some suggestive features. 
We postulated an equation that balances cyclotron-resonant and non-cyclotron-resonant 
dissipation effects of kinetic origin with steady spectral transfer into the dissipation range 
due to MHD scale cascade processes. This formal structure evidently has been able to 



- 11 - 



account for some of the observed properties of the distribution of inertial range cross helicity 
and dissipation range magnetic helicities. We note that several important approximations 
are implicit in our treatment. For example, we do not account in any way for the energy in 
the velocity field, E v in the dissipation range. This might be an acceptable approximation in 
either of two cases: if E v and its dissipation rate are much smaller than E b and its cascade 
rate S; or if a proportionality or approximate equality exists between E v and E^ in the 
dissipation range. In the absence of a better theoretical guidance, as well as plasma data 
at the requisite frequencies, we prefer the latter explanation at present. We can suppose, 
for example, that the "Alfven effect" attempts to enforce near equipartition of velocity and 
magnetic fields, or that the phenomenological dissipation rates we assumed in fact include 
some contributions that are mediated by couplings to the velocity field, which would be 
expected to be heavily damped at scales near the thermal proton gyroradius. 

Motivated by the typical values of magnetic helicity in the dissipation range (Fig. 1), 
and assuming that all fluctuations propagate in one direction, we estimated near equality 
of cyclotron resonant contributions represented by 7 r and other dissipative mechanisms 
represented by 70 . Using the inertial range cross helicity to estimate the relative likelihood 
of propagation direction produced a refined estimate j r ~ 270. 

In this development we have been forced, due to limitations of spacecraft 
instrumentation, to use the inertial range cross-helicity to compute a proxy for the 
propagation direction in the dissipation range. The most notable limitation of this 
substitution derives from the possibility that preferential dissipation may lead to different 
cross-helicity values in the dissipation range, although we are not aware of any observational 
evidence for this. Indeed, the connection between a c and direction of propagation may be 
complicated in the dissipation range by various kinetic wave modes such as whistlers. On 
the other hand lower frequency observations of a c (see Matthaeus & Goldstein 1982) often 
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indicate that a single direction of propagation is dominant over several orders of magnitude 
of scale, which would tend to support our extrapolation into the higher dissipation range 
frequencies. In any case, the correlation evident in Fig. 1 appears encouraging with regard 
to use of this proxy. 

The present results provide some preliminary insights into the structure of the 
interplanetary dissipation range, but additional work needs to be done to better understand 
the physics of the kinetic dissipation mechanisms represented by 70 and j r . For example, 
we expect on general grounds, that Landau and non-resonant processes should make a 
contribution to dissipation of three-dimensional, MHD turbulent fluctuations, but an 
acceptable large amplitude theory of such processes is not yet developed as far as we are 
aware. Similarly, resonant dissipation, generally evaluated by linear Vlasov theory, requires 
improvement for the same reasons. In addition, linear theory makes no prediction about 
damping of purely transverse "two-dimensional" turbulence, which appears to be favored 
by MHD in the presence of a moderately strong mean magnetic field (Matthaeus et al. 
1994). In this regard one would expect that MHD turbulence would be accompanied by a 
turbulent induced electric field E = — v x b that would produce stochastic acceleration of 
suprathermal particles, and associated damping of the fluctuations. Further developments 
in kinetic theory are required to describe dissipation that is nonlinear, anisotropic and 
driven by an MHD cascade. 
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Fig. 1. — Scatter plot, for 33 Wind data intervals, of the normalized cross helicity in the 
inertial range, a c , vs. the normalized magnetic helicity in the dissipation range, a m . Triangles 
are intervals with outward directed mean magnetic field, and bullets have inwards mean fields. 
The dashed line corresponds to the best-fit line through the origin, a c = — 1.90<7 m . 
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